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The long-time behavior of the velocity autocorrelation function in a classical two- 
dimensional electric conduction system is studied by the molecular dynamics simulation. 
In equilibrium, the effect of coexistence of many-body interactions and a random poten- 
(f) , tial is investigated. A crossover from a positive tail proportional to t~ , to a negative tail 

proportional to — 1~ 2 is observed as the strength of the random potential increases. In 
nonequilibrium, the positive tail is enhanced whereas the negative tail appears at earlier 
times as an electric field increases. 



§1. Introduction 



The long-time tail was first discovered by B. J. Alder and T. E. Wainwright 1 ' for 

^ ■ the velocity autocorrelation function (VACF) in a hard-core fluid system. This tail 

is positive and proportional to t~ ' 2 (hereafter called the fluid- type tail). 2 ' 7 > Here, 

d is the dimension of the system. As a result the self-diffusion coefficient D in this 

system is logarithmically divergent in terms of the system size L for d = 2. 

Another system which has a long-time tail is the Lorentz model, which describes 
the motion of a single particle in a disordered system. In this system the tail of 
the VACF is negative and proportional to — £~("+ 2 )/ 2 (called the Lorentz-type tail 
hereafter). 8 ' 13 ^ Although D is not divergent even for d = 2, this tail might cause 
an algebraically L-dependent term in D. 
J^T ■ In the hard-core fluid only a many-body interaction exists and in the Lorentz 

[^. \ model only a random potential exists. Hence a question occurs: How does the long- 

time tail change when a many -body interaction and a random potential coexist? A 
. ■ typical system which has both a many-body interaction and a random potential 

is an electric conduction system. In this system electron-electron (e-e) interaction 
corresponds to a many-body interaction and electron- impurity (e-i) interaction to a 
^^ . random potential. 

In studies of the long-time tails in the hard-core fluid and the Lorentz model 
the equilibrium states have mainly been investigated. Regarding electric conduction 
systems, another interesting question arises: How does the long-time tail change 
under a nonequilibrium condition? 

In this paper we report the results of the molecular dynamics (MD) simulations 
on a model of a two-dimensional electric conduction system to answer these two 
questions. 14 ) 
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§2. Model and its physical meanings 

We use a model of electric conduction proposed in Ref. 15). We here explain its 
physical meanings in more detail. 

The model is a two-dimensional classical system, the size of which is L x x L y . 
In the system are three sorts of particles, which we call electrons, phonons and 
impurities. 16 ) An external electric field E is applied in the x-direction, which acts 
only on electrons. The boundary condition in the x-direction is periodic, and that in 
the y-direction is a potential wall for electrons and a thermal wall with temperature 
T for phonons. Moreover we assume that short-range interactions are present among 
all these particles. The charge of an electron is denoted by e. 1G > The mass, radius, 
number density of the electrons (phonons) are denoted by m e (m p ), R e (R p ) and 
n e (n p ), respectively. The impurities are immobile and play the role of a random 
potential. The radius and number density of the impurities are denoted by Ri and 
rii, respectively. The configuration of the impurities is given by an almost uniform 
distribution except for the restriction that the distance between any pair of the 
impurities is larger than 2(R e + Ri). The initial positions of the electrons and 
phonos are randomly arranged not to overlap with the other particles, and their 
initial velocities are given by the Maxwell distribution with temperature T. 

We can control the frequencies of the e-e and e-i collisions by changing n e and 
m. Thus we can investigate independently the effects of the many-body interactions 
(e-e, e-p and p-p interactions) and that of the random potential (e-i interaction) 
on the long-time behavior of the VACF. It should be also noted that the situation 
where n, = corresponds to the hard-disk fluid and that the situation without e-e 
interaction corresponds to the (non-overlapping) Lorentz model. 

A typical system corresponding to this model is a two-dimensional electron sys- 
tem in a doped semiconductor at room temperature. Since the Fermi energy is 
smaller than room temperature in such a system, electrons can be treated as classi- 
cal particles. Furthermore, we can treat them as being confined in a two-dimensional 
plane because the temperature is lower than the exciting energy between the ground 
and the second subbands 17 ) of an electron. 

It is well-known that in a uniform solid (such as our model) , the motion of the 
electron system can be separated into the collective oscillations and the individual 
motions, and that the long-range effects of the Coulomb force are adequately incor- 
porated in the collective oscillations. 18 ^ Moreover the collective oscillations do not 
play the central role in the electric conduction. [It is worth mentioning here that for 
non-uniform systems, such as a conductor connected to electron reservoirs, 19 )' 20 ) the 
long-range effects should be treated more carefully as discussed in Refs. 19), 20).] For 
these reasons, we can treat the electrons as individual particles interacting through 
a short-range force when discussing transport properties of conductors. We can es- 
timate the range (the screening length) of the effective interaction among electrons 
by the Debye length, (A;BT'/47rn e e 2 ) 1 ' 2 , in such a classical system. The potential 
produced by an impurity is also screened and its effective potential range is also 
estimated by the Debye length. Therefore the interaction ranges of e-e and e-i in- 
teractions are comparable. 
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In real solids phonons are the oscillation modes of the crystal lattice in the 
conductors. Therefore the total number of phonons does not conserve. However, 
because almost all the possible modes of phonos are excited in semiconductor at room 
temperature, the number density of phonons is so high that the non-conservation of 
the phonon number would be irrelevant. Moreover, the energy-momentum dispersion 
relations of phonons are complicated in real solids. This would be also irrelevant, 
however, when discussing general nonequilibrium properties, which are independent 
of the details of the materials, of electric conduction. In this study, therefore, it 
is sufficient to model phonons as classical particles 21 ) whose number conserves and 
whose mass is constant, which corresponds to parabolic dispersion relation. Note 
that if one wants to reproduce the T-dependence of conductivity, such as the famous 
linear-T dependence, then the number of phonons should be varied as a function of 
T. In this paper, however, we are not interested in the T-dependence of conductivity. 

This model contains what we believe to be essential elements of electric con- 
duction in the context of nonequilibrium statistical physics. These elements are 
the following: (i) A driving force which induces electric current: electric field, (ii) 
Careers which transfer heat to outside the conductor: phonons. In real physical sys- 
tems (i.e., in experiments), a conductor is surrounded by a large insulating material, 
which works as the heat bath. The energy supplied from an external electric field 
to electrons is dissipated as the Joule heat transferring into the heat bath through 
the walls of the conductor to keep the system in a steady state. The heat flow 
across the walls of the conductor is mediated not by electrons but by phonons, 
while the heat flow in the conductor is mediated by both of them. Therefore the 
electron-phonon interactions and heat contact of phonons at the walls are essential 
to realize a nonequilibrium steady state, (iii) Objects which break the microscopic 
translational invariance in the bulk region: impurities. Although the conductor is 
macroscopically uniform, the microscopic translation invariance is violated by the 
imperfections of the conductor (impurities, defects and so on). This defines the rest 
frame of the electrons at equilibrium, and thereby eliminates any possible anoma- 
lies which may arise from the translational invariance. (iv) Many-body interactions: 
short-range interactions among all particles, (v) A nonequilibrium steady state is 
uniquely determined by a small number of macroscopic parameters: T and E. 

We perform the time-step-driven MD on this model. We set R e , m e , k-gT and 
e to unity, and take L x = 2000, L y = 125, R p = 1, Ri = 0.5 and m p = 1 in 
these units. 23 ) The number densities of the electrons and phonons are fixed to 
n e = 0.04 and n p = 0.002. Since n e is taken larger than that in our previous 
study, 14 ) the long-time tails are more clearly observed. In this study we employ the 
Hertzian interaction [7™ 1 = y(max{0, djk}) 5 ' 2 as the interaction potential between 
any pair of the particles. Here, Y is a constant (fixed to 4000 in the simulation), 
and djk = Rj + Rk — \ r j — r k\ is the overlap of the potential ranges of the j-th and 
A;-th particles {Rj is the radius of the j-th particle and equals R e , R p or Ri. rj is 
the position of the j-th particle). 23 ) 

We previously confirmed that this model shows good properties of steady states. 15 ) 
For example, as shown in Fig. [TJ we observe linear response of v x to E near equi- 
librium and nonlinear response far from equilibrium, where v x is the component 
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Fig. 1. Average velocities of an electron for (a) m = 0.0004 and (b) m = 0.04, plotted against E. 

of electron's velocity in the x-direction (parallel to E). We also observed that 
the Kramers-Kronig relation of the complex admittance as well as the fluctuation- 
dissipation theorem 22 ^ hold at all frequencies in the equilibrium state. 

§3. Results 



We calculate the VACF, 



C{t) 



[Mt) - <^))M0) - («,») 



((v x - K» 2 > 
of an electron in steady states under various conditions. 



(3-1) 



3.1. Equilibrium case 

First we show the results of the VACF in the equilibrium states. We calculate 
the VACFs with n e and n p fixed, and rij varied from to 0.04. In this case, we 
can investigate the effect of the random potential introduced into the fluid system. 
Figure [2] displays the results. We observe power-law tails of the VACFs at longer 
times. When rii = the tail is positive and the exponent is approximately — l, 24 ) 
which is the result for the fluid system. As rij increases, the fluid-type tail becomes 
weaker, and eventually disappears. Instead, we observe another tail which is negative 
and whose exponent is approximately —2. That is, the long-time behavior of the 
VACF is a crossover from the fluid-type tail to the Lorentz-type tail as n, becomes 
larger. 

According to phenomenological theories, 2 -*' 4 ' 7 '' 12 '' 13 ' the origin of the long-time 
tail in a fluid is the mode coupling of the momentum field and the density field (be- 
cause of the conservation of the total momentum and total number of the particles) 
while that in the Lorentz model is the coupling of the density field and a static mode 
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Fig. 2. A double logarithmic plot of the absolute value of the VACF in the equilibrium states for 
various values of m. The open and closed symbols correspond to positive and negative values 
of the VACF, respectively. The dotted and dashed lines are reference lines proportional to t _1 
and t~ 2 , respectively. 



generated by the configuration of the impurities. From these theories, we have a 
qualitative explanation of the crossover from the fluid-type tail to the Lorentz-type 
tail. When m is small, the momentum conservation is approximately valid and we 
can regard the system as a fluid. Then the momentum field can contribute to the 
hydrodynamic mode and we observe the fluid-type tail. As rii increases, the vio- 
lation of the momentum conservation becomes larger and simultaneously the static 
mode by the impurity configuration becomes more relevant. Then the fluid- type tail 
disappears and the Lorentz-type tail appears. 

A necessary condition to treat the electron system as a fluid is that the mean 
number N e of the electrons in an area among the impurities is sufficiently large. Be- 
cause the mean area per impurity is l/n», N e is estimated as N e = n e /ni. Therefore 
we can use the ratio of the electron and impurity density as a criterion for the long- 
time behavior of the VACF. From Fig. [2] the crossover value (n e /nj) c is evaluated 
as 2 < (n e /r&j) c < 10 [the fluid-type tail appears when n e /nj > (n e /ni) c and the 
Lorentz-type tail appears when n e /ni < (n e /ni) c ]. 



3.2. Nonequilibrium case 

Next, we show the results of the VACF when electric current is flowing. 



We 



calculate the VACFs for n e /ni ^> (n e /rii) c and n e /ni <C (n e /rii) c , with varying E. 
In Fig. [3l we show the long-time behavior of the VACFs of the steady states including 
the states in the nonlinear response regimes as well as those in equilibrium and in 
the linear response regimes. When n^ = 0.0004 [n e /rii ^$> (n e /nj) c ], we observe that 
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Fig. 3. A double logarithmic plot of the absolute value of the VACF in the steady states for various 
values of E. The open and closed symbols correspond to positive and negative values of the 
VACF, respectively. The dotted and dashed lines are reference lines proportional to t _1 and t~ 2 , 
respectively, (a) m — 0.0004 \n e /rii S> (n e /rii) c ; a corresponding E-(v x ) plot is in Fig. Ufa)]. 
(b) m = 0.04 [n e /ni <C (n e /n,:) c ; a corresponding E-{v x ) plot is in Fig. [Tfb)]. 



the amplitude of the fluid-type tail gets enhanced as E increases. Because C(t) is 
normalized by ((y x — (i^)) 2 ) to be C(0) = 1, this enhancement is not simply due 
to the rise of the kinetic temperature of the electrons. This result implies that a 
"temperature" is different among different time scales in nonequilibrium states. For 
the system with raj = 0.04 [ra e /n, <C (n e /ni) c ], the Lorentz-type tail appears at 
earlier times as E becomes larger. One of the reasons for the early appearance of 
the Lorentz-type tail would be that the duration time of the e-i collision becomes 
shorter as E increases. 



§4. Concluding Remarks 



In summary, we have investigated the long-time behavior of the VACF in a sys- 
tem where many-body interactions and a random potential coexist. In equilibrium, 
we have observed a crossover from the fluid-type tail to the Lorentz-type tail as 
the impurity density increases. We have interpreted that this crossover occurs be- 
cause the system can not be regarded as a fluid when the impurity density becomes 
large. The ratio of the electron and impurity densities is a criterion quantity for this 
crossover. In nonequilibrium, as an electric field increases we have observed that 
the fluid- type tail is enhanced for n e /rii S> (ra e /nj) c and that the Lorentz-type tail 
appears at earlier times for n e /rii <C (n e /rii) c . 

Finally, we list some issues related to this study. 

(1) The long-time tails in equilibrium states induce system-size dependence of D. 
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In an electric conduction system D can be translated into the electrical conductivity 
a in the linear response regime by the Einstein relation. Thus some system-size 
dependence of a should be observed in experiments. To our knowledge, however, 
no experiment has been reported which observed a system-size dependence of a for 
uniform two-dimensional electron systems at room temperature. 25 ' This might be 
because the amplitude of the system-size dependent term in a is too small to detect 
in comparison with the system-size independent term, or because phonon scattering 
would introduce a system-size independent cutoff. 

(2) In Sec. 3.1 we have presented a possible scenario of a mode-coupling theory 
for the crossover from the fluid-type tail to the Lorentz-type tail in equilibrium. 
This should be explicitly shown by combining the mode-coupling theories for the 
fluid 2 '' 4 ' 7 ' and for the Lorentz model. 12 '' 13 ' 

(3) This crossover behavior would also be supported by a kinetic theory with 
a two-parameter (n e and rij) expansion, which combines the kinetic theory for the 
hard-core fluid (n e -expansion) 3 ' and that for the Lorentz gas (n^-expansion). 8 ' 

(4) A crossover might occur from the fluid-type tail to the Lorentz-type tail, 
even in a system with lower density of impurities [that is, n e /ni > (n e /nj) c ], at 
longer times when the effect of the violation of the momentum conservation becomes 
relevant. If this is true, D is convergent 26 ' in two-dimensional systems except when 
the impurity density vanishes. To demonstrate this, a larger-scale simulation with 
higher accuracy is necessary. This issue should also be studied by a mode-coupling 
theory and a kinetic theory. 

(5) In Sec. 3.2 we have observed an enhancement of the fluid-type tail in a 
nonequilibrium steady state. This might be explained by a mode-coupling theory 
similar to the one recently developed in sheared fluids. 27 ' 
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